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ABSTRACT 


This  report,  consisting  o(  Parts  I,  II  and  III,  is  concerned  with  three 
different  aspects  of  the  response  of  arches  and  donees  under  dynamic  loads. 

In  Part  I  che  accuracy  of  an  approximate  design  method  for  arches  sub¬ 
jected  to  dynamic  loads  is  evaluated  by  comparing  the  predictions  of  this  method 
with  the  exact  solutions.  Primary  emphasis  is  placed  on  the  effects  of  loads 
v^ich  are  uniformly  distributed  around  the  arch. 

In  Part  II  the  response  of  circular  elastic  arches  under  a  moving  pressure 
pulse  is  investigated  by  the  modal  method  of  analysis.  Various  combinations  of 
natural  modes  are  considered,  and  the  combination  of  the  smallest  number  of 
modes  which  satisfactorily  approximates  the  exact  solution  is  determined.  It 
is  concluded  that  good  approximation  to  the  true  response  can  be  obtained  by 
considering  the  contributions  of  the  first  antisymmetrical  and  the  first  two  sym¬ 
metrical  natural  modes  of  vibration. 

Part  III  presents  a  derivation  of  an  approximate  theory  for  the  dynamic 
response  of  spherical  shells  loaded  unsymmctrically  by  time-varying  pressures. 
Non- linear  effects  are  considered  so  that  the  resulting  equations  reflect  the 
buckling  tendencies  associated  with  large-deftection  behavior. 
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SVALUATION  OF  A  ISSIGN  MBTECD  FOH  ASCHSS  SUBJECTED  TO  DYNAMIC  LOADS 


1.  Cb.ieet  and  Scope 

Til?  purpose  of  the  stuiy  dessribcd  in  this  part  of  the  report 
was  to  evaluate  the  accuracy  of  the  design  oethod  presented  as  Reference  1* 
for  arches  subjected  to  the  pressures  associated  with  an  atomic  explosion. 

Ihe  approach  used  consists  of  evaluating  the  response  of  several  arches  by 
this  aiethod  and  coaparing  the  results  vith  those  obtained  on  the  basis  of 
the  method  reported  in  Reference  2,  the  latter  method^  the  actual  arch 
is  replaced  by  ft  discrete  framework  consisting  of  a  series  of  rigid  bars; 
flexible  Joints^  and  concentrated  point  masses.  For  the  sake  of  brevity^ 
this'method  will  be  referred  to  as  the  "exact"  method,  and  the  method  of 
Ref.  1  as  the  approximate  method. 

The  solutions  presented  are  for  two-hinged,  circular,  elastic 
arches  of  uniform  cross  section.  IRie  majority  of  the  solutions  are  for  u 
uniform  all-aroxind  pressure  vith  a  tli^-vlse  variation  represented  by  a 
triangle  with  on  initial  peak.  Some  resiilto  are  also  presented  for  a  tri¬ 
angular  pressure  pulse  moving  across  tlie  arch. 

Inasmuch  as  Reference  1  is  not  z^adlly  available,  a  brief  de¬ 
scription  of  the  approximate  method  is  given  in  the  following  section.  The 
"exact"  method  is  described  in  oetail  in  Reference  2* 

Review  of  Approximate  Method 

In  this  method  the  actual  loading  on  the  arch  Is  separated  Into 
a  symmetrical  component  and  an  antlsymmetrical  component,  as  shown  In  Fig.  1. 
The  intensity  of  the  pressure  for  the  syanetrical  component,  p^.  Is  con¬ 
sidered  to  be  constant  around  the  arch.  Siatll^ly,  for  tne  untlsysmeUical 
component,  the  pressure  intensities,  for  ti»-*  vindvarc  ond  leeward  sides 
are  considered  to  be  equal  but  of  opposite  si^i. 

♦  Refereiiue  flgui'e  and  equation  numbers  refer  to  those  in  each  part  of  the 
report. 


Ihe  rcsj'onee  unaer  the  c>ainetrlcel  component  of  loading  is  evalu¬ 
ated  considering  tvo  modes  of  defomation: 

(a)  Ibe  "uniform  ccspnjssion**  node,  ^ich  is  considered  to 
produce  a  unlfona  compression  vithout  an>'  bending,  and 

(b)  Ibe  "bending-cotipression*  mode  in  which  the  arch  bends  in¬ 
ward  at  the  crown  and  outverd  near  the  haunches.  It  is 
assumed  that  the  bending  moment  at  the  crown  is  due  to 
this  mode  only. 

nie  intensity  of  the  load  associated  with  the  imifom  compression  mode  is 
taken  as  the  intensity  of  the  symmetrical  load  component,  whereas  tlie  in¬ 
tensity  of  the  load  associated  with  the  compression-bending  mode  is  taken 
as  one-third  of  that  for  the  unlfora  c<wpre88ion  mode. 

Ihe  response  due  to  the  entisyaaetrical  component  of  loading  is 
evaluated  by  considering  the  so-called  "deflection  node"  which  is  antl- 
syometrical  about  the  crown.  Thus  the  computation  of  the  response  in  each 
mode  involves  the  analysis  of  a  sln47l)-degree-of-freedom  system. 

The  natural  periods  of  the  arch  iji  the  various  modes  of  deforma¬ 
tion  are  determined  as  follows:  For  a  circtilar  arch  with  uniform  cross 
section,  the  natural  period  of  the  unifom  compression  mode,  is  taken 
as 


where  a  denotes  the  mass  per  unit  of  length  of  the  arch,  R  the  radius  of 
the  arch,  A  the  cross  sectional  area,  and  E  the  i^dulus  of  elasticity  of 
the  material  of  which  tlie  arch  is  caponed. 

The  period  of  the  bending-compressior.  mode  for  no  thrust  in  the 
arch,  is  taken  equal  to  the  period  of  e  slnnly  supported  beva  h'lving 
a  span  equal  to  one-third  the  a.,  /.cn^th  of  the  arch,  i.e. 
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synbol  9^  denotes  the  angle  of  opening  of  the  arch>  as  shovn  in  ?lg.  1. 
5he  prlae  on  is  used  to  indicate  that  the  effect  of  the  thrust  is  not 
considered  in  this  expression. 

The  corresponding  period  of  the  deflection  code,  is  taken 
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and  is  a  correction  factor  accounting  ror  the  fact  that  the  arch  does 
not  act  entirely  os  a  beau.  The  follovlng  expressions  are  given  for 
For  a  hinged  archj 

2  .  - 

(6a> 


and  for  a  fixed  arch 


where 
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The  "actual"  perlodSi  and  ^d'  wnlch  Jncox’porcte  the  effect 
of  the  thrust,  are  deteralned  approxlsatcly  hv  .'•'»''.tip3ying  the  values  of 
T*  and  Ti  by  the  fact*.T 
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>*ere  p  is  the  intensity  of  the  "model  pressure"  on  the  ai’ch  av  the  in¬ 
stant  that  the  corresponding  node!  displac^aent  is  saxlxaios  and  is  the 
critical  static  huck3ing  pressure.  For  a  arch,  p^^  is  given  by  the 

equation 


With  the  intensity  of  the  loading  and  the  value  of  the  natural 
period  or  vibration  for  each  mode  of  defoxnation  determined,  the  response 
of  the  arch  in  each  of  these  redes  is  ev^uated  "iry  analyzing  the  arch  as  a 
system  with  a  single  degree  of  freedom. 

For  an  arch  subjected  to  the  ejasiietrical  component  of  loading 
shown  in  Pig.  1,  the  axial  force  associated  with  the  uniform  compression 
mode  is  determined  from  the  equation 

H.  (A.F.)^p^R,  (9) 

where  p^K  is  approximately  the  force  produced  imder  static  conditions,  and 
(A.F.  is  the  amplification  factor  for  R.  ^Is  factor  depends  on  the  shape 
of  the  applied  pressure  and  the  ratio  of  ISaC  duration  of  the  pulse  to  the 
natural  period  cf  vibration  of  the  particular  mode  of  deformation  considered. 
The  maximum  .bending  moment  associated  with  bending-compression  mode 
occurs  at  the  quarter  points  of  the  arch*  and  ma;'  be  expressed  as 


M.  (A.F.)„H, 


(10) 
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viierc  the  static  sooentj  is  taken  as 


■  p^Rr  . 
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1  -  (3?) 
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Ihe  first  factor  on  the  rig^t  side  of  this  expression  represents  the  oaxi* 
Bi«  Boaent  in  a  slsply  supported  beaa  bavtrg  a  span  length  as  defined 
by  Sq.  (3)>  and  subjected  to  a  unifom  static  pressure  of  Intensity  p^/3* 
the  syabol  r  represents  the  radius  of  gyration  of  the  cross  section  of  the 
arch.  ?^e  second  factor  accounts  for  the  effect  of  arch  cur/ature.  The 
aaplificatlon  factor  (A.P.)|^  oust  be  deterained  on  the  basis  of  Uie  period 
of  the  coapresslon-beiuiing  mode,  T*^. 

For  the  antlsyanetrical  component  of  loading  shown  in  Fig.  X,  the 
bending  noaect  is  deter&ined  froa  £q.  (lO).  However,  the  aspllflcatlon 
factor  must  be  based  on  the  period  or  the  deflection  sodc  T^,  and  the 
static  Doaent  must  be  detexvlned  for  a  pressure  p^  and  a  length  as  given 
by  Eq.  (3)<  For  a  hinged  arch,  the  e:pxe68ion  for  becomes 


32 
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(12) 


Ihe  effect  of  the  axial  thrust  is  taken  into  account  in  figuring 
the  bending  resistance  of  the  arch.  For  the  compression  bct.dlng  mode,  the 
bending  re^iXstance  is  taken  as  the  product  of  the  corrcspondLng  resistance 
for  no  thrust  multiplied  by  the  factor 

1  -  p/pjp 

where  Is  the  buckling  load  corresponding  to  a  length  -  R  For 

the  deflection  mode,  the  reduction  factor  is  V*o-  as 

1  -  p/p.. 
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■*ere  is  as  previovaly  defined.  5he  Inverse  of  these  factors  can  also 
he  interpreted  as  aagaiflcEtlon  factors  to  be  applied  to  Bj.  (lo). 

'•  gc.r.~-lts  of  Ccatarlsons  for  Arches  tbder  Ustfors  Pressure 

3.1  General.  Ihe  lesponse  carves  presented  la  She  reaaladtr  of 
this  part  of  the  report  are  for  tvo-binged  circolar  elastic  arches  of  imi- 
foxn  cress  section  subjected  to  a  unifom  all-around  pressure  pulse.  3he 
pressure-tlae  relationship  Is  represented  hy  an  Initially  peaked  triangle, 
as  ehowi  In  the  upper  left  comer  of  Fig.  2a.  Ihe  peak  intensity  of  the 
pressure  Is  denoted  by  p^  and  the  duration  of  the  pulse  by  t^.  n,e  values 
of  these  quantities  are  specified  in  tenns  of  the  dlaenslonless  ratios 
^o^^cr  ^d^^o'  ^cr  t^Ptresents  the  critical  buckling  pressure  cor¬ 

responding  to  an  antlsymtetrlcal  node  of  deformation,  and  represents  the 
breathing  period  of  vibration  of  a  ccetplete  ring  having  the  same  radius  and 
cross  section  as  the  orche 

It  has  been  shovn  (Ref.  2)  that,  for  the  conditions  considered, 
the  arlal  force  In  the  arch  Is  fsirly  uniform  along  the  arch  and  that  the 
maximum  moment  occurs  at  the  crown.  In  view  of  this,  only  the  moment  and 
the  axial  force  at  the  crown  will  be  considered. 

Since  the  arch  la  uniformly  loaded,  la  the  application  of  the 
approximate  me'uiod,  onlj-  the  uniform  compression  mode  and  the  bendlng- 
coapresolon  mode  need  be  coosidered.  Ihe  former  mode  yle’.i'  the  axial  force 
in  the  arch,  and  the  latter  gives  the  mooient  at  the  crown.  For  tl.e  tri¬ 
angular  pulse  Investigated,  the  aoplificatlon  factor,  A.F.,  for  an  clastic 
system  that  is  Initially  at  rest  la  given  hy  the  following  equations: 

For  t  <  t^, 

A.?.  '  1  -  ^  -  cos  2«  I  4  ^  u,  e.  I  (Ij) 
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For  t  > 

I  t  -  * 

A.?.  «  'SjA®  +  sin  (2«  — a)  (lU) 

’-Ixere 

A  •  COB  2*  ^  ^  sto  2S  (llib) 

d 

B  «  -  (1  -  cos  2*  y)  +  =ln2*^  (I’w) 

d 

a  >  Un"^  I  (ihd) 

O 

In  the  above  ejq)ression8^  the  quantity  T  denotes  the  natural  period  of  the 
particular  node  of  defoxnation  considered. 

5.2  Comparison  of  Typical  Response  Curves.  Li  Fig.  2  are  shovn 
the  tiae*historles  of  the  axial  forxe  and  the  bending  moment  at  the  crovn  of 
a  tvo'hlnged  arch  with  the  folloving  dimensions: 

Bise-to-Span  Hatlo  «  f/l.^  «  0.£0  (corresponds  to  9*^  »  87-21°) 
Slenderness  Ratio  «  L^/r  ■  100  (corresponds  to  B/r  «  72-5) 

Tilt'  intensity  of  the  peak  presstire  and  the  duration  of  the  pulse 
are  as  follows: 

Vl“cr  ■  ^  V’b  ■  = 

The  dashed  '^urve,  reproduced  from  Reference  2,  may  be  considered  to  represent 
the  exact  solution.  The  solid  curve  was  obtained  by  the  approxizrtate  method- 
From  Fig.  2a  it  ceji  be  seen  that^  except  for  a  slight  phase  dlfler' 
ence  ^icb  is  of  no  importance  from  a  design  point  of  view,  the  approximate 
solution  for  the  axial  force  is  in  very  good  ev-^r'^eni  with  the  "exact"  solu- 
ticn.  Included  in  this  figure  is  also  the  soluticL  obtained  by  xc>lacing  the 
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valuc:  ol  the  nattiral  period  of  the  usifors  compression  acdc  by  the  value  of 
the  "extenslonal  mode*  of  vlbr&tion  of  the  arch.  For  the  particular  arch 
Investlffated,  the  extenslonal  node  corz^pox«ls  to  the  second  s^raetrlcal 
node  of  vibration  (see  Fig.  4.2b  of  Refer^ce  2^.  She  use  of  this  period 
does  not  alter  the  solution  significantly. 

In  Fig.  2b^  althou^  the  peak  values  of  the  bending  ncnent  de* 
ternined  ^  the  tvo  methods  are  in  reasonable  agreement^  the  periods  of 
the  oscillations  are  significantly  different.  She  indications  are  that  the 
approxiaate  method  does  not  acco\mt  properly  for  the  behavior  ol  the  systea 
In  bending. 

In  order  to  gain  better  insight  Into  the  true  behavior  of  the 
system^  the  moment  at  the  crovn  vas  evaluated  ty  application  of  the  iscKial 
nethod  of  analysis*  In  ^Is  analysis,  the  effect  of  the  axial  thrust  on 
bending  vas  neglected.  In  Fig.  3  are  given  the  modal  contributions  for  the 
first  three  symnetrlcal  modes.  In  addition,  the  sm  of  these  contributions 
Is  compared  with  the  exact  solution.  It  can  be  seen  that  the  solution  based 
on  the  first  three  modes  Is  In  good  agreement  vltn  the  exact  solution.  She 
Indications  are  that  all  three  contributions,  particularly  the  one  asso¬ 
ciated  with  the  second  or  extenslonal  zsodc,  are  quite  important. 

3*3  Effect  of  Pressure  Intensity  Parameter.  In  Figs.  4  and  S  are 
shown  time  histories  of  the  axial  force  and  the  be.'vding  moment  at  the  crown 
of  the  arch  considered  before,  except  that  the  Intensity  of  tnc  peak  prescurc 
has  the  values  of  Pq/p^j.  "0.5  and  2.0.  She  other  parameters  of 

the  problem  are  the  same  es  those  Indicated  in  Fig.  2.  She  "exact"  solu¬ 
tions  presented  in  these  and  all  subsequent  flgitrec  are  ieproduccd>  from 
Reference  2.  In  Figs.  4a  and  5a,  the  solid  curves  are  Identical  wlt!i  Use 
corresponding  curve  In  Fig.  2a,  since  in  the  **pp.*oxlmate  method  the  axial 


force  is  considered  to  be  inieperdent  cf  the  ratio  it  is  of  interest 

to  note  that  the  best  aereenent  between  the  apprcxixaate  and  exact  solutions 
for  axial  force  is  obtained  for  Fig.  idii^  corresposds  to  the  ssallest 
value  of  p  /p  considered.  In  general,  the  obser/ations  th^t  can  be  mace 
about  the  cur/es  in  Figs,  k  and  5  ere  similar  to  those  made  about  Fig.  2. 

^.4  Effect  of  Pulse  duration  Parameter.  ‘She  arch  considered  in 
the  preceding  section  was  also  analyzed  for  values  of  t^/T^  in  tlie  range 
between  0.25  and  4.0.  In  Fig.  6  are  ijiven  the  time-histories  of  the  axial 
force  and  bending  moment  at  the  crown  for  a  value  of  ®  0*75#  and  in 

Fig.  7  are  given  spectrum  curves  for  the  complete  range  of  t^/T^  values  con¬ 
sidered.  These  curves  express  tlie  absolute  maximum  value  of  the  axial  force 
and  bending  sKxaent  at  the  crown  as  a  function  of  the  parameter  t^/T'^* 
the  exact  and  the  approximate  solutions  are  indicated.  It  can  be  seen  that 
the  Diaximum  axial  force  predicted  by  the  dppxt>ximate  method  is  in  very  good 
agreement  with  that  obtained  by  the  "exact"  method.  On  the  other  hand, 
there  are  significant  differences  between  the  values  of  tlie  maximum  moments 
predicted  by  the  two  methods. 

3.5  Hffect  of  Arch  Dimensions.  Comparative  solutions  were  ob¬ 
tained  for  a  nimibcr  of  arches  with  different  dimensions*  Two  groups  of 
problems  were  considered: 

(a)  Arches  with  a  rise-span  ratio,  of  0.2,  end  slenderness 

ratios,  l*^/r,  of  50,  100  and  200. 

(b)  Arches  with  a  value  of  L^/r  •  100  and  rise-span  ratios  of 
0.1,  0.2  and  0.5.  The  latter  ratios  correspond  to  values  of  eoual  to 
45. 24*^,  97*21®  and  I80®,  respectively.  The  Iod-I  parameterc  ore  taken  us 

Po/p„  ’  i  “'5  V^O  ■ 
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results  are  eisuaried  in  Pigs,  d  tbrou^  13*  can  ts  seen 
that^  in  all  csses^  the  peak  value  of  the  axial  force  predicte^l  by  the 
approxlmave'oetnod  is  in  excellent  agreessent  vlth  the  exact  value.  On  the 
other  hand^  both  the  tise  histories  and  the  peak  values  of  the  bending 
a^oments  corresponding  to  the  tvo  methods  of  solution  differ  significantly 
from  one  another,  peek  values  of  the  b^^nding  moments  are  compared  in 
Pigs.  10  and  13 . 

J*.  Comparisons  for  Arches  Subjected  to  a  Moving  Pressure  Pulse 

Xn  Pig.  ll)  are  shovn  the  results  obtained  for  an  arch  subjected  to 
a  moving  pressure  pxdse.  Oie  pressure-tine  relationship  vas  represented  by 
an  initially  peeked  triangle.  Ihe  problem  parameters  were: 

f/L^  .  O.a  L^r  .  100 

V^cr  “  ‘A  ■  ^ 

‘Zhe  quantity  t^  represents  the  transit  tlme^  l.e.  tlic  time  required  for  the 
front  of  the  pulse  to  move  across  the  arch.  !nie  remaining  symbols  have  the 
same  meaning  as  before. 

Xn  the  approximate  solution,  the  peek  intensity  of  the  symmetrical 
component  of  loading  (See  Fig.  1}  was  taken  as 


Pc - Po  '  ^^5) 

1.5  *  O.jf 

and  the  corresponding  intensity  of  the  antisymmetricel  or  defies tlon-cwdc 
component  vao  taken  as 


Pd 


0.5  + 0.5 -f 
1.5  +  0.5  i 


Po 


(16) 


The  pressure-time  relationship  and  the  duration  of  t.iC6e  component  Juadlngs 
were  considered  to  be  the  same  as  those  of  tne  actual  pulse. 


?3roTi  the  res'ilts  presented  ia  Pig.  15  it  can  he  seen  that  the 
approximate  solution  for  the  axial  for^e  is  not  in  as  good  agreement  vith 
the  exact  solution  as  for  the  cases  eonsiucreo  previously.  It  may  be  vorth 
nothing^  however^  that  in  this  e^se  the  stress  correspondl’zg  to  the  maximum 
social  force  is  relatively  smell  in  ccotparlson  to  that  resulting  from  the 
maximum  bending  moment.  Ihe  maximum  hecdi.x,j  I'cxaent  occurs  at  the  5/^  pointy 
and  its  peah  value  is  predicted  vith  reasonable  accuracy^  as  can  be  seen 
free  Fig,  l4c. 

In  Fig.  15  are  given  response  curves  for  moment  at  the  l/k  and 
5A  points  of  the  arch  considered  before,  except  that  the  dxiration  of  the 
pulse  is  represented  by  a  value  of  ■  2.  It  can  be  seen  that  the  agree¬ 

ment  between  the  solutions  obtained  ty  the  two  methods  is  fairly  good,  par¬ 
ticularly  ^dien  tlie  curves  corresporrflng  to  the  absolute  maximum  effect  are 
considered. 

5. 

Ihe  results  of  the  comparat..ve  stti^les  presented  can  be  summarized 
as  follows: 

(1)  For  arches  subjected  to  a  uniform  radial  pressure,  the  approxi¬ 
mate  method  can  predict  vith  excellent  accuracy  the  magnitude  of  the  maximum 
axial  force  in  the  arch.  However,  both  the  time  history  and  the  magnitude 
of  the  maximum  xooaent  determined  ty  this  method  may  be  sigr.lficantly  dlflerent 
from  the  exact  results.  The  results  of  a  modal  analysis  that  has  been  mode 
suggest  that  the  computation  of  the  bending  laoment  at  the  crown  requires 
consideration  of  the  first  three  syametricel  modes  of  vibration. 

(2}  On  the  basis  of  Uie  two  solutions  obtained  for  arches  sub¬ 
jected  to  a  triangular  moving  pressure,  it  that  the  value  of  the 

maximum  moment  predicted  by  the  approximate  method  Is  in  reasont.ole  agree¬ 
ment  vith  the  exact  value.  TV.e  axial  forces  are  not  pivdicted  as  accurately 
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in  this  case  as  for  a  syacetricai  loa^;  however^  this  saay  DOt  be  very 
significant  since  the  stress  resulting  froa  the  maxlnum  axial  force  is 
usually  small  in  cccQ)arison  to  that  due  vo  the  aaxiaua  noaent. 

(^)  It  is  believed  that,  for  elastic  arches,  a  more  rational 
approximate  method  of  analysis  can  be  developed  by  considering  the  con¬ 
tributions  of  the  first  tvo  or  three  fyiroptrlcal  nodes  of  vibration  and 
of  the  first  antisyometrlcal  mode. 
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DYNAMIC  ANALYSIS  OP  CIRCULAR  ELASTIC  ARCHES  BY  HOBAL  MEIHOD 


i,  Ob.lect  and  Scope 

The  general  purpose  or  this  study  was  to  investigate  -lie  dynaalc 
response  of  elastic  arches  by  application  of  the  raodal  method  of  analysis. 
Specifically,  it  vas  desired  to  evali^te  the  rolative  contribution  of  the 
various  natural  modes  to  the  total  response,  and  to  establish  a  basis  for 
the  development  of  a  simplified  mett»>d  of  anal^'sis. 

Numerical  solutions  have  been  obtained  for  c  range  of  the  para¬ 
meters  involved,  and  the  complete  solutions  incorporating  the  contributions 
of  all  the  natural  s^des  of  vibration  verc  compared  vith  corresponding 
approximate  solutions  obtained  by  considering  a  limited  number  of  modes.  An 
attempt  vas  made  to  determine  the  combination  of  tlie  smallest  number  of  modes 
which  best  approximates  the  complete  solution.  All  solutions  are  for  tvo- 
hlnged,  circular,  elastic  arches  subjected  to  a  moving  triangular  pressure 
pulse.  In  the  analysis,  the  buckling  tendencies  of  the  arch  are  neglected; 
accoi^ingly,  the  response  *s  proportional  to  Uie  intensity  of  the  applied 
pressure. 

A  digital  computer  program  vas  developed  for  the  numerical  computa¬ 
tion  of  the  response.  This  program  has  teen  developed  specifically  for  an 
arch  subjected  to  a  moving  pressure  pulse  of  triangular  shape.  The  replace¬ 
ment  system  is  conside}«d  to  have  ten  bars  of  equal  horizontal  projections. 

A  variable  mode-counter  set  in  the  program  makes  it  possible  to  consider  any 
desired  nicaber  of  modes. 

In  addition  to  the  numerical  results,  brief  descriptions  of  the 
substitute  structure  used  in  the  method  of  analytes  are  Included  in  tills 


report. 
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2.  Description  of  the  Model 

A  general  view  of  the  physical  uodeX  used  to  approximate  the 
actu»al  arch  is  shovn  in  Fig.  1.  model  la  similar  to  that  presented  in 
Ref.  1  j  with  the  exception  that  the  continuous  arch  is  repl^.c'»d  by  ten 
bars  of  equal  horizontal  projectioi^  instead  of  ten  bars  of  equal  length. 
With  this  arrangement  the  computation  of  the  Instantaneous  distribution  of 
the  load  over  the  arch  is  simplified.  Ihe  bars  are  considered  to  be  laass- 
less  and  rigid  in  both  bending  and  axial  deformation;  the  displacements  aitd 
angle  charges  are  defined  at  the  Joints. 

Tbe  lumped  mass  at  any  Joint  J  is  given  by  the  equation 


J  -AfX 


irtiere  is  the  length  of  the  Jth  bar  and  p;  is  the  mass  per  unit  length  of 
the  continuous  arch. 

Similarly  the  stiffness  of  the  flexible  Joint  at  J  is  derived  from 
the  elastic  properties  of  the  actual  arch  between  mid-points  of  consecutive 
panels  and  is  given  by  the  equation 


J 

where  is  the  stiffness  of  Joint  J,  E  is  the  elastic  modulus,  and  Z  is  the 
moment  of  Inertia  of  the  cross  section  of  the  arch. 

bars  of  the  analogous  framework  are  numbered  consecutively 
starting  with  J  ■  0  at  the  left  hinge  and  tensinatlng  with  J  •  10  at  the  right 


*  *'Re6ponsc  of  Arches  under  Dynamic  loads,"  by  H.T.  Epplnk  and  A.S.  Vel'T.tsoi., 
Uhiverslty  of  Ulinois  Project,  AFSW  -  TR-^-55/  Air  Force  Special  Weapons 
Center^  Klrtland  Air  Pbrcc  Base,  Rev  Mexico,  I>*ccmber  5950. 
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hinge^  as  shovn  In  Fig*  !•  Bars  arc  numbered  In  the  same  order  Trom  1  to  10, 
with  the  bar  connecting  joints  j-1  and  j  designated  as  bar  j.  The  total 
central  opening  angle  of  tiie  arch  Js  denoted  by  9^.  ^e  deformed  configura¬ 
tion  of  the  arch  axis  is  defined  with  respect  to  the  undefoiucd  shape  iii 
terns  of  the  radial  and  tangential  cooponents  of  the  joint  displacements. 

For  joint  j  these  are  designated  as  w,  and  v,,  respcctiveJy.  Tlte  displace- 

V  « 

ment  v  is  considered  as  positive  radially  inward,  and  v.  positive  clockwise, 
J  J 

as  indicated  in  Fig.  1.  Ihe  rotations  of  the  bars  are  considered  positive 
in  a  clockwise  sense. 

3,  Method  of  Analysis 

Ihe  method  of  analysis  is  based  on  the  modal  analysis  technique 
as  applied  to  the  substitute  structure,  ^e  data  required  for  such  an  analy¬ 
sis  include  the  natural  frequencies  and  modes  of  vibration  of  the  model. 

Since  the  distribution  of  the  pressure  along  the  arch  varies  with  time,  the 
participation  factors  for  the  various  modes  are  time-dependent  quantities. 

Ihe  analysis  presented  here  is  restricted  to  elastic  range,  end  buckling 
tendencies  are  ignored.  In  other  words,  results  are  valid  only  for  sjaall 
values  of  where  is  the  buckling  pressure  of  the  arch. 

Ihe  computation  of  the  response  involves  the  following  steps: 
i.  Determination  of  the  natural  frequencies  and  modes  of  the  sub¬ 
stitute  siiuCUue. 

11.  Determination  of  the  load  distribution  across  the  arch  at  any 

time. 

ill.  Computation  of  the  modal  participation  factors  corresponding 
to  a  discrete  system  of  time  Intervals,  the  time  increment  chosen  beuig  small 
enou^  to  yield  solutions  which  ai'e  reaeonfd>V‘  o'tcuiatxi.  Such  c  computation 
Involves  the  derivation  of  the  governing  partlcjjatlon  factor  x.^..aticn  for 


each  mode.  An  apprcxiuale  marching  f.pe  Af  procedure  vas  adopted  for  tlie 
solution  of  the  par'Jiclpation  factor  equation.  In  this  procedure,  the 
Instantaneous  participation  factors  and  their  time  derivatives  for  each 
mode  are  obtained  by  use  of  the  appropriate  quantities  at  ihc  previous  ♦^iae 
instant. 

iv.  Superposition  of  the  modes  at  any  time  Instant  considering 
the  respective  participation  factors. 

In  this  report  the  details  of  the  procedures  are  omitted.  All 
solutions  were  obtained  on  the  IIXIAC,  the  digital  compviter  of  the  Jnlver- 
sity  of  Illinois. 

1*.  Numerical  Solutions 

4.1  Problem  Parameters.  Ihe  properties  of  the  arch  arc  specified 
in  terns  of  the  folloving  parameters: 

•  the  total  opening  angle  of  the  arch. 

L^/r  ■  the  ratio  of  the  opan  of  the  arch  to  the  radius  of  gyra¬ 
tion  of  the  cross-section  of  the  arch,  referred  to  as  the 
slenderness  ratio. 

Ihe  load  parameters  are  defined  with  respect  to  the  fundamental 
breathir.g  period,  T^,  of  a  complete  ring  of  the  same  radius  and  thickness. 
This  period  is  given  by  the  equation 


Ihc  lopd  parameters  aic 

^  luratlon  of  the  blast  pulse 

"  Breathing  period  of  the  complete  ring 

\  Time  of  travel  of  the  peak  rre*suse  front  over  th«*  arc):  span 
*  Breathing  period  of  the  complete  zing 


4.2  Frobleas  Considered.  <iicerical  solutions  vere  obtained  for 


the  foHoving  values  of  the  paraceters: 


’0 

■  60° , 

lJt  -  100, 

‘A  - 

-  2  and  5 

’’0 

■  &pf 

50, 

‘A  * 

-  2  ana  5 

***0 

-  90°, 

L  /r  -  100, 

0  ' 

‘A  - 

=  3  and  7.5 

^0 

-  90®, 

yr-  50, 

.  1.5, 

yi^  .  3  end  7.5 

In  all  cases  the  time  Interval  of  Inteeration  vas  taken  as  O.05  T^. 

O^e  response  of  the  arch  vas  evaluated  for  the  following  condl' 

blona: 

(i)  Considering  all  eighteen  nodes  of  vibration 
(11)  considering  the  first  antls^Boetrlcal  mode  and  the  first  and 
second  sysaietrical  nodes.  Ihe  latter  solution  will  be  referred  to  as  the 
throe-fflode  solution.  In  addition,  for  several  of  the  problocs,  solutions 
were  obtained  considering  only  two  oodcs  of  vibration.  The  ccsbinatlons 
considered  included  the  first  syioaetrlcal  aj.d  the  first  antis>Ts:etrical  nodes, 
or  Uic  second  symaetrical  and  the  first  antlsyscctrlcol  nodes.  Finally, 
some  solutions  were  obtained  by  considering  the  first  two  syisctrlcal  and  the 
first  two  antlsyxBcctrical  modes.  In  this  report  only  a  few  of  these  solu¬ 
tions  are  presented. 

4.5  Comparison  of  Solutions  Obtained  by  Two  Different  MeUiodc.  As 
a  check  on  the  accuracy  of  the  computer  progres  used  to  obtain  the  num^^lcol 
data  presented  herein.  In  Figs.  ^  and  U  the  time  histories  of  the  dloplocc- 
aents,  DOuents,  and  thrusts  detemined  by  the  sixteen  node  solution  ore  com¬ 
pared  with  those  determined  by  application  of  the  aettiod  described  In  the 
iefeiciKe  listed  on  V.  2.2*  These  results  are  «  /  an  arch  with  9.  »  60^  and 
h^/r  «  100,  T^ie  first  four  circular  natural  frequencies  and  t*»*»  corresponding 
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laodes  of  vibr&tios  of  this  si'ch  are  sc^vn  in  Tig.  2.  It  can  he  seen  that 
there  is  fairly  ^od  agreement  between  the  two  solutions^  the  Ezaxisun 
difierenee  for  the  peak  values  of  displ^ements,  ooments,  and  thrusts 
being  less  than  ^  percent.  !Ihe  phase  difference  between  the  two  solutions 
is  due  to  the  sli^tly  different  procedures  employed  in  the  two  methods  tc 
concentrate  the  external  pressure  at  the  node  points. 

4.4  Presentation  and  Dlscnsstoa  of  Results.  For  eaca  of  the 
problems  enumerated  abeve,  the  solutions  were  studied  to  determine  the 
combinaticn  of  nodes  which  best  approximated  the  sixteen  mode  solution. 

Since  the  arch  model  considered  has  eighteen  degrees  of  freedom;  the 
eighteen  mode  solution  represents  the  exact  solution,  "in  general;  it  was 
fc'ond  that  the  solution  based  on  the  first  antisymetrlcal  mode  end  the 
first  two  symmetrical  modes  provided  a  very  good  approximation  to  the  true 
response. 

In  Figs.  5  through  14  are  given  the  time  histories  of  the  response 
curves  for  displacements,  moments;  end  cxlal  forces  at  selected  points  of 
the  arch  for  several  of  the  problnss  considered.  Both  the  three  mode  solu¬ 
tion  and  the  ei^teen  mode  solution  are  presented,  it  can  be  seen  that  the 
two  sets  of  solutions  are  consistently  in  good  agreement.  !Zhe  best  agreement 
is  obtained  for  displacements  and  axi^  forces.  Ihe  moment  curves  are  in¬ 
fluenced  to  a  greater  extent  by  the  hl^er  modes  than  eithCii  .ne  deflection 
or  tiie  axial  force  c\irves. 

For  displacements;  reasonable  agreement  was  also  obtained  by  con¬ 
sidering  only  the  first  synstrical  and  the  first  antisyumetrical  modes. 
However,  the  results  for  moment  and  axial  forces  were  generally  unsatiolactory. 
^he  solutions  based  on  the  first  aatisyaactriev  .code  and  the  second  sya- 


raetrical  sode  vere  generally  no  better  ^en  the  tvo  laode  solution  r^x'erred 
to  shove.  For  displacements  and  xaoaents,  the  greatest  contribution  to  the 
response  ^as  provided  by  the  first  antlsyaoetrlcal  node. 

5.  Stggy^ary 

Ihe  results  of  this  study  indicate  that  a  satisfactory  method  of 
analysis  for  arches  txnder  moving  pressures  can  be  developed  by  considering 
only  the  first  antlsyirtetrlcal  and  the  first  tvo  symmetrical  natural  modes 
of  vibration. 
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AN  APPROXIMATE  THEORY  FOR  TER  LARGE  DYKAJIIC  DEFIPCTIONS 
\  OF  UNSIMMBTRICALLY  LOADED  THIN  SPHERICAL  SHELLS 

Introduction 

Presented  herein  is  a  general  tl^ry  for  the  large  deflections,  but 
nnall  strains,  of  spherical  shells  leaded  uncyriAtrically  by  time-varying, 
external  pressures.  Three  eqiuatlons  in  three  displacement  unknowns  are 
developed  by  a  variational  procedure  ^ploying  Hamilton's  Principle  of 
Least  Action. 

The  strain  energy  function  used  in  the  formulation  of  the  total 
energy  of  the  shell  is  developed  tader  the  sasumption  tha*.  ibc-  state  of 
stress  is  approximately  plane;  that  is,  the  effect  of  transverse  shear  stress 
and  of  transverse  normal  stress,  acting  on  surfaces  parallel  to  the  middle 
surface,  may  be  neglected.  Thia  assrsnption  yields  the  usual  Hooke's  Lav 
relations  used  in  ordinary  plate  theory. 

The  lav  of  deformation  of  a  fiver  of  the  shell  is  governed  by  the 
Kirchoff-Love  Assumptions  vbich  sta*«  that  {a)  points  lying  on  a  normal  to 
the  xmdefomed  middle  surface  remain  on  the  same  normal  to  the  deformed 
middle  surface,  and  (b)  the  dIeplaCemeits  in  the  direction  of  the  normal  to 
the  middle  surface  are  approximately  equal  for  all  points  on  the  some  normal. 

These  assumptions  restrict  the  foU.ovirig  analysis  u  thin  shells 
<  1^) .  However,  the  equations  of  motion  include  the  effects  of  rotatory 
inertia. 

In  the  following  development,  the  external  pressure  is  assumed  to 
be  caused  by  air*  blast,  and  is  always  directed  normal  to  the  outer  surface 
of  the  shell. 

The  significant  tontributloae  of  this  stuJy  include  the  i'>rtcntion 
of  all  the  important  nonlinear  terms  in  the  strain-displacement  relations, 
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aad  a  rigorous  calcxilation  of  the  potential  energy  of  the  external  pressure 
to  sake  possible  the  consideration  of  "dynamic  irstability." 


Exact  Strain-Displacement  KelalionH 

The  complete  stiain-disp^ccaent  relations  for  a  spherical  deme 

(2> 

have  been  derived  earlier,'  *  and  oxe  svmarized  here  for  convenient  reference. 
The  position  of  any  point  on  the  sheli  in  its  undefomed  state  is  defined  by 
the  orthogonal  curvilinear  coordinates  (9,V,s).  .^The  displnccmcmi  vector 
(u,Y,v)  is  chosen  such  that  u  and  v  are  tangent  to  tnc  Q  and  coordinate 
lines,  vhile  v  is  directed  outvard  along  the  surface  noxroal  (see  Fig.  l), 

in  the  following  relations,  subscripts  on  u,  v,  v  indicate  partial 
derivatives,  while  subscripts  on  €  and  7  Indicate  the  direction  of  the  strain 
component.  The  radial  distance  to  a  point  of  the  shell,  a  +  s,  is  denoted  by 
r. 


■M' 


*  cot  9  +  v) 


(:) 


’’87.  -  <“2 

r  r' 

r 

[(«j)(V'')  *  (''zHVg)  +  («j,)(Wg-U)j 

’'92  • 

> 

--) 

r' 

r  sin  0 

-  V 

cot  8)  +  +  u  cot  9  +  u)  + 

Assmaptlon  about  Displacement  Variation  Through  Thickness 

'The  usual  assmption  in  plates  and  shells  chat  normals  remain 

normal;  straight  and  inextenslcnal  in  the  deformed  state  is  used  here.  This 
assumption  about  middle  surface  normals  leads  to  the  result  that  the  displace¬ 
ment  vector  (u,V;w)  varies  linearly  with  z  through  the  shell  thicJaiess.  Thus 
ve  have: 

u  (6,<P,z)  u  Uq  (e,<p,o)  +  2  ig  (6,^,0) 

V  (9,<P,2)  -  Vg  (S,II>,0)  +  26^  (2) 

w  (9,'P,2)  »  Wg  <8,9,0) 

where  (UqjV^jV^)  is  the  diCiilacemeuv  vector  of  a  point  on  the  middle  surface 
at  s  »  0.  is  the  rotation  of  the  noiaol  in  the  meridional  plane.  5’^  is  the 
rotation  of  this  normal  In  a  plane  normal  to  the  meridional  plane.  If  the 
rotations  of  the  normal  are  snail;  w(z)  will  not  differ  markedly  freo  v(z»0); 
and  so  V  is  assuned  to  be  independent  of  t  ond  &. 

The  rotations  of  the  normal  are  determined  from  the  Kirchofr-L'ive 
Assumptiorut.  For  the  normal  to  remain  straight  during  deformation^  no  shear 
defonoatlon  zaay  take  place.  This  requires  that 

’'02“’'92'® 

SuliSvitutlon  of  Equations  (1)  and  (2)  Into  Eq.  \  yields  two  olnultoncous 
equations  relating  9,  5  oa-.d  their  derivatives. 


<Jf  (a-i.v+«y*5|rg)  +  8  (Vg+zSg)  -  u  - 


*  (^^-6  -  -  V  cot  e  -  z  S  cot  9) 


{5a) 


V  2  6  V 

*  ®  +  u  cot  6  +  t  fr  cot  9  +  w)  .  V  -  -jf- 


Since  ijr  and  6  are  independent  of  z,  their  values  are  required  to 
satisfy  the  auxiliary  conditions:  *  tif^  +  66^  ■  0.  An  approxinate 

solu^-ion  to  Eqi  (3e)  is 


a+w+trt:ot  &  +  cp 

sin  9 


(5b) 


finally,  to  facilitate  the  derivation  of  the  equations  of  notion, 
we  approxiaate  Eq.  (5b)  as 


t 


a 


sin  $ 


Approxinnte  Straln-Dlsplaceaent  Relations 

Fpcu  the  three-dimensional  expression  for  Hooke’s  Lav,  the  require¬ 
ment  that  0,  •  0  results  in 
z 

't  ■  ^  's  ^  V  j 

vMch  gives  h'  ^ 

Thus  it  5  s  necessary  only  to  evaluate  €^,  and  7^  directly.  Substitution 
of  Sqs.  (2)  and  (4)  into  (l)  results  in  (retaining  only  nonlinear  tcxns  ir.  v 
and  its  derivatives): 


-^•5- 


-  (-2  -  ) 
0  'a  a+z  a  a+z' 


if, 


(6) 


— ^  1  (v  -  v„)^  +  (a+z)2(-i)2  cot  0)'' ! 

2(a«)2  0  ^  ej.  <?  /  j 


%-(■ 


u  cot  Q 


a  sin  9  a-f-z  a-fz 


2(a+z)‘ 


■  [(V  -  I  Wg 


'  sin^  9 


) 


(7) 


Jk 


e  V  cot  e  2z 


'09  'a  sia  0  a 


a>z  a  ain  0  a+z  a  sin  0 


+  -tr — r-  (v.  cot  0  -  v„,)^ 

O  try  j 


22_!£i!l!) 


cot  9  w, 


9q> 


2(a+2)^  L  a*^sia  9 

V’  cot  9  w, 


3  -  I  ''00^^  -  sll^)  (8) 

W  sin  9  a  sin  9'^''  a  2.  a  '^9 


The  term  l/a^z  taay  he  expanded  as: 


1  i  Z  ^  Z  X 

a+z  **  a  ^  *  a  „2  ' 


/  i  \  2  1  /»  /X  z  I  z  » 

W  “:2(l-2a-^5-2-  •••) 


and  the  principal  strains  uoy  be  expressed  in  the  fom: 

Cg  .  (A  +  Bz  +  Cz-  f  ...) 

»  (D  +  Fz  *  Hz?  +  ...) 


(9) 

(10) 

(11) 


(J  *  Kz-i  Li  +  ...) 


Fo21owir.g  the  suggestion  of  Laagbaar,^^^  the  expansions  for  the  strains  are 
cut  off  aflw*r  the  .Unear  tem  In  a.  Suhstitution  of  Eqs*  (9)  and  (lO)  into 
Bqs.  (6),  (7)  end  (8)  then  yields  (rctainlz:^  only  derivatives  of  w): 


A-f  (u^  +  w) 

B  -  -  ^  (v  +  w^g) 


1  ,  <D®  t 

F  =  .  —  (  w  +  Wg  cot  e  »  ■ 


^  •  7  (iir  * ''0  - *■ 


(12) 


Stress  Resultants 

In  the  case  of  plane  stress,  the  elastic  stress-strain  relations 
assuae  the  fora: 


”9  ■  ‘‘P’  '  ’“P  "  **<?  *  '’  ‘9’’  V  ■  ’’W 


UJ) 


Substitution  of  Eqs.  (11)  and  (IJ)  into  the  foliovi?.^  stress  resultant 

foraulas  results  in: 

h 

«0  ■  r  "9  "  f)  0'  ■  t 

Oh  1-V  *•  J 

*  2 

h 

f  (i  ■*“  f)  [(I>*-'A)  ■*■  (FfvB)  j  (15) 

’  2 
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T 


%  ■  'V  ■  ’e,  ^  *  fe  («)]  (16) 

5  (A  +  vD)j  (17) 


Mg  .  r  2  ZOg  (1  +  |)  dz  -  — [(b  +  VF) 
■'ll  ®  lafi-v"^)  L 


la(i-v^) 

2 

-  3  - 

'V  ■  .  i  (i>.vA)  j 

2 

-  5 

“a?  ■  ‘V  ■  fl  ^^69  + 1)  *  ai^r:7  [''♦ !('’)] 


(18) 


(19) 


h 
’  2 

Strain  Energy 

?or  fitooll  strains,  the  elastic  strain  energy  of  deformation  is 
given  by  the  volume  Integral: 

“•  I  /^V  [«a«a  *  v?  *  ® 


this  can  be  expressed  in  terms  of  strains  only  by  use  of  Bqs  (15). 

h 
K 


U  ' 


^  2 


(2i) 


(a'fz)  sin  9  d0d9dz 
Replacing  the  strains  by  Eqs.  (U)  and  performing  the  integration,  wc  hav^: 


U  «  +  Ug  + 


vhere: 


■'2x  r  9  I 


2(l-v‘^)  -'O  '^0  ' 

v.hJ„2  /'2*  /■?  I 


(£8) 

']  ait.  0  d0dp 

(23) 

1^)  K®j 

1  sin  0  d0d? 

m 
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Eh^ 


|a®+D®+2iAD  +  (■^)  J®1 

„  C  sin  ? 

+  tej^AB+DF+v{.\F+BD)+(^)JKj| 


d0cH> 


(25) 


is  the  meabrane  strain  energy  due  to  stretching  of  the  middle  surface. 

U.5  is  the  usual  expression  for  the  flexural  acruln  energy  used  in  most  theories 
of  thin  shells.  is  a  correction  to  the  luending  energy  due  to  the  Inclusion 
of  the  teir  (a+z)^  instead  of  in  Eq.  (21). 


Kinetic  Energy 

The  kinetic  energy  is  given  by  the  volume  integral 
r2j(  pe 
■'0 


-irrrh  +  z+^)^  +  (v^  +  sin  0  d0d<pdz 

2  ^  (26) 


Potential  Energy  of  External  Pressure 

During  deformation;  the  elastic  stiain  energy  changes.  AlsO;  the 
external  presrure  does  vork,  and  the  corresponding  change  of  energy  oust  be 
included  in  the  expression  for  the  total  energy. 

The  external  pressure  is  considered  to  be  caused  by  air  blast*  As 
such  it  is  directed  normal  to  the  outer  surface  of  the  shell.  The  potential 
energy  of  the  external  pressure  is  then  given  as: 

R  ■  -  JJf’i  ^<^'’'1 

vhere  is  the  change  of  an  Infinitesimal  element  cf  volume.  Since  the 
gectaetry  of  the  shell  is  specified  in  polar  coordinates  It  Is  convenient 

to  use  these  coordinates  in  speciiying  the  volu:*  .ngc.  An  infinitesimal 
element  of  shell  volume  is  expressed  as: 

dV  .  P®  ain  dSg  dp 
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vhere  p  is  the  distance  free  the  center  of  usdefomed  cur/ature  to  a  point  on 
the  outer  surface.  Calling  the  distance  to  a  point  on  the  defomed  outer 
surface*  and  the  distance  to  a  p^iat  on  the  undefomed  outer  surface,  the 
potential  energ/  of  the  external  pressure  becemes: 


1  r^-'r  r-j 
5  Jo  Jo  r  ■> 


sin  Sj,  d9jj  asPjj  -  py  sin  0  dSdipl 


(28) 


vhere : 


pQ  -  -\/(a  + 1  +  +  (u*)®  +  (v*)^ 


(29) 


and  the  angles  locating  the  defomed  point  are  approximately : 


Oj)  "  9  +  [  (- 


•  9  » 


l-Sa  +  5  f 


...1 


-)  - 


(30) 


-  - /  5\  u - / 

g  r  V*)  sin  $  u*  cos  0  (a  +  ^  +  w^^lsin  0  +  u»  cos  0 


...] 


vhere  the  starred  quantities  indicate  displacements  on  the  outer  surface, 

2  »  Equations  (20)  indicate  that  0^  and  9^^  do  not  depend  on  the  squares 
of  the  displacements.  Thus,  vc  retain  only  the  linear  teraa  these  expres¬ 
sions,  and  have: 


d0n 


(19„  »  49  + 


[nfh:] 

tinTr 


fio 


(51) 


L(a  +  2  *"  ^ 


These  relations  are  approximate^  ond  based  on  the  assxnption  that  0|^  is  a 
function  of  0  only,  and  9^  is  a  function  of  9  only.  With  the  relation  tliat 


-i.iO- 


?in  0  »  sluj  B  +  - - —  I  o  sla  sl  coaf - - 

°  t  a  +  ^  *  v»J  L  ^  +  J 


-> 


(J2) 


+  cot  0  sin  ('  '  “  ■“■)  I 

a  +  I  +  vx  J 

the  potential  energy  express  on  heccaes: 

^  "  5  y'  y'  ^  ^  ^  ^  ^ 


05) 


vhere:  x  ■  • 


+  I  +  V*  a  V  I 


(a  +  ^  +  v*)sin  0  +  M*  cos  6  ^  ^ 


Finally  tha  integrand  aay  be  expanded  by  ualng  the  following  approxisationo: 

J 

j^(a  +  I  H  v»)^  t  (u»)^  +  (v»)^J  I  (a  +  I  <•  v»)^  +  |(a  +  5  +  w»)j^(u»)^i(v*)®j 


.  U#  *  ^  1  y 

COS  X  ■  cos  C  ^  )  ■  1  •  2  ' 


u*  k2  ^  1  /  u*  \2 

-ij-— )  -1-5  (—ij) 


sin  X  ■  sin  (* 


'  h  *•  h 

n+^  +  v*  a+-  +  v  ^'*’2 


«B*  ♦  V/sin  9 


X  ■\  V  ^  U^*  ♦  V 

a  t  r;  la  +  ^joln  B  a  +  - 


la  t  il®  eln  a 
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Sul)stitution  of  these  approxiaate  relaticas  into  Eq.  (35)  results  la  the 
foliovirv  expression  for  the  lioteatlal  eaeigy  (w*  •  v): 


”  ^fo  Jo  ^  *  f  *  -  (a  +  |)  sin  9  dSd?  (}L) 


h  2  ,  2  ,  9 

"  '  2  ^  2  1 


vhere:  h.  =  (a  +  5  +  v)  +5  (a»)  ^  ^  (y») 


,  u*  sot  9  1  /  n*  , 

!'o  •  J-  * - S - o  ( - h) 


a  +  | 


(55) 

(56) 


'V,*' 
_2 _ 


.  ,  V  -  ^  V _ 

"5"^  a. I  (a.|)2aln9 


_2_ 


(37) 


Equations  of  Motion 

The  principle  of  virtual  vork,  when  extended  to  a  dynoalcs  problea 

yields  the  useful  relation 
t, 

(er  -  5y  -  sn)  dt  -  0  (i8) 

This  is  often  referred  to  as  ”HaffliItoii*s  Principle.”  For  a  conservative 
system^  Eq.  (38)  becoaes 

0  /  (T  -  y  -  n)  dt  •  0  (39) 

j''0 


L 


Equation  (39)  iaplles  that  enong  all  actions  that  will  carry  a  conservative 
systea  froa  a  given  configuration  to  a  given  configuration  in  a  given 
time  interval  (t^it^),  that  which  actually  provides  a  stationary  value 

to  the  integral.  The  slsnlflcance  of  Eq.  (39)  i*'  that  the  Sulur-;xigraage 
equations  for  the  integral  are  the  differential  equations  of  aotlon  for  the 
systea  under  ecnslderatlon. 
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Substitutloa  of  energy  ftoecions  from  Eqs.  (23),  (zU),  (25),  (26), 
and  (3^)  -Uto  Eq.  (33)  yields  the  follovlag  variational  equation: 

_  _ -  >2  ,  .21. 

2. 


5  ri  ^  r  “  2  2-1 

2  Jj,  *'  J  '] 


ofir^T' 

Jq  •''O 


-g-  i^A®  ^  C®  +  2;*D  t  <-^^)  j2  j 


Eh\2  (-^2  ^  „2 

2lt(l-v2) 


+  2vBF  + 


(^>K2j 


7sln  e  dfldWt  «  0 


2ll(l-»®) 


^a2  +  f-  +  2vAD  +  (^)  j2^ 

[  t  <«  -^AB  +  DF  +  V  (AF  +  BD)  +  (■^)  J!^ 


5  *  I  *  ”)('‘iV3^  ■  ^  “  ■*■  |)^ 


/ 


(iO) 


vhich  becomes 


ff7\ 

Jfc  ^0  ''0  \ 


s  fix  t  C0(u2  ,  v2)  ,  (i  ,  n)(.,2)  , 

I  sin  0 

L  Si 


-  g^^  --  I  a2  .  ..2  ^  ^  .1-A  ^1.0 

9  A  +  D  +  2vAD  +  .*»  i  sin  9 

2(l-v‘^)  I  2  J 


l; 


EhV  j 

+  F^  +  2vBF  + 

(^)  K' 

■]  sin 

2M1-v‘^)  I 

L 

J 

FhJ 

'"ra  2 

>tK  +  +  2vAD 

+  (^-) 

2'.(l-v®) 

1 

'  2  ' 

+  ^  r  AB  +  DF  ■• 

V  (A"  V 

BD)  + 

fhi 

t  vlCUjhjkj)  -  (a 

fin  0 

/d^dWtnO 


n 


slu  9 


j  w 
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12a‘ 


2  * 


Designating  the  integrand  oC  Eq.  (Ul)  Ir/  the  sitsbol  “f,"  the 
Euler'Lagrange  Equations  assuae  the  fom: 

bt  d  /5f  1  d  jdt  \  ()  rt 

'  -55  ■  -59  ^  ° 


dr  d  /df  \  d  /df  %  d  /df  \  * 

^  -  ■3a  -  -35  '3r>  -  3t  <Sr>  ‘  ° 


(te) 

(W) 


df  d/df\  o  j  VI  \  V  ivi  \  .  V /oi  »  <j  /  oi  \ 

33  -  ^  <3^>  - -53  (33;)  *  ^  (3;^)  " (3;^) 


d  ,df  , 


d  ,df  , 


/df  _ 


df  ■ 


m 


dt' 


89 


9t 


Performing  the  operations  indicated  in  Eas.  (42),  {^5)^  and  (44)  results  in 
three  nonlinear  partial  differential  equations  la  three  unknowns  (u,v,w). 
Many  of  the  terms  in  these  eqi^tions  nay  i>e  replaced  by  their  syabcHc 
representation  as  stress  resultants*  These  thrie  equations^  in  terms  of 
stress  resultants  are  presented  belov. 


35 


0  "0 


Mj 


^  h  ^  V  *  ®  ^”0  -  ^  ■  V. 


r  1  r  3o  3  3  .^'1 

■  [“tt  ^  5)  -  T-J  -  UTe  [  ^  -  35  ^  '33;^ 

^  (-“V  ^  »69>  "  3lb  I  ^  V  ^  2  cot  s  .  M  J 


m 


(W) 


"9tt  1  ^  1  fbq  _  d  ^ 

If  Gin  0J  *  sin  0  ^  ^ 
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^  Ul.  0  (H,  *  «.«,)]  *  ^  (M,  -  a=K,)] 

+  2  [^'e?  "■  ““atJ  ■  ®  [“  ■"  V  ■"  ®  '  V  j 

5  r'*©  ^  ^.>  .  _«  \  •  *P  '•'  *  nW.  ^  -  f.os  0  ("M-  +  anlJ».)l 

+  "Se  [ — I —  ^‘‘e  ■*■  "  ^Tsin  e  "0V  V'  ^ 

- 1  [  ^  "  2^  <“<»? "  “'W  "  ^ 


o  ph^  sin  0 


■»  1  Vtt  \ 

■tt<^  "■  fe)  ^  ^  “  *  “ett  ^ 

h 


.-i5<''0tt<=°‘®"''e8u^;^> 


-  ft  - 1?  ^  <%>  ^  fa  ^  fa  <%;>  ^  ^  ] 


uh?re:  Q  ®  (a  +  ^  + 

In  the  limit,  vhcre  a— ^0,  and  vhcrc  the  nonlinearity  of  Zv  (**?)  i® 
neglected,  these  equations  ^ssme  the  fora: 

h  i<-\  *  «ei  "  »ib  i  ^”^8?  "  "  'O'-  ®  y^o  ■  V  ■  '”8  ■  "y  ] 


(■‘7) 

(W) 


a  p  ha**  u,*. 


('•s) 

+  K„J  »  Pha'’  ',v,  (50) 


I  *  V' "  ^  h  "  V  * "  ®  I  “”89 "  H '•« 

^  ^  ,  cPt  8  -  cPt 

■  9  ■*■  o  .2  2la  0  <30  dp  ^ 

^0*^  sin  0  op 


(5l) 
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ate  left-hand  sides  of  Eqs.  (ItJ),  (50),  and  (jlj  are  in  the  sane  fora  as 
those  dc''ived  by  by  static  cqullib»*im  considerations. 

Equations  (^5)  through  (h7),  vhen  coupled  vith  stress  resultant 
foitrulas  yield  three  equations  in  three  unknowns.  These  equations  Include 
nonlinear,  large-displaceaent  effects,  rotatory  inertia  effects,  and  the 
effects  of  change  of  geooietry  on  the  dynaolc  eqimibriica  of  the  systes. 

this  change  of  geometry  effect  enters  in  two  vays.  It  enters  all 
three  equations  through  the  energy  of  the  external  pressure  on  a  shell 
surface  whose  area  changes  with  tine.  It  also  enters  the  eqmtion  for  radial 
equilibrius,  2q.  (47),  due  to  the  large  deformation  effect  which  alters  the 
line  of  action  of  the  stress  resultants.  This  latter  effect  is  observed  only 
in  Eq.  (47)  because  only  nonlinear  terms  in  w  and  its  derivatives  have  been 
retained  in  Fqs.  (12). 

The  nonlinear  equations  ore  too  lengthy  to  express  here  in  displace¬ 
ment  form,  however,  the  linearized  Bqs.  (49)  through  (51)  are  presented  here: 


v„cot  e 


Sin  B 


*  •<'«  ( 


1  r  «  cot  0  T 
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Equations  (52)  through  (Sh)  agree  exactly  uith  Eqs.  85(a-c)  of  Ref.  (3)  except 
for  the  omission  of  a  term  in  Eq.  85(a)  of  Flujge's  text.  These  equations  arc 
adequate  to  use  in  the  coapucation  of  the  natural  frequencies  of  the  spherical 
shell.  To  find  the  buoklirg  prsosureo,  houcYsr,  the  nonlinear  teins  of 
Eqs.  (45)  through  (47)  must  he  added  to  these  equations.  Also,  the  rotatory 
inertia  toms  con  no  longer  he  emitted. 


-5.17- 


REFEB&ICES 


(1)  H.  L.  ianghaar,  "A  Strain-Energy  Expression  for  Thin  Elastic  Shells,’* 
Jonmil  of  Applied  MechmUcs,  p.  I85,  JOne,  1949, 

(2)  J.  B.  Ihiberg,  R.  L.  Jennings,  et  al.,  "Analysis  and  Design  of  Domes, 
Arches,  and  Shells,"  Vol.  IX  of  Interim  Reijort  on  Contract  No, 

AF  29(601) '?591,  PjroJect  No.  lOSo^  Dept,  of  Civil  Engineering, 
University  of  Illinois,  p.  C-6,  Noveaccs,  i960. 

(3)  W.  Flugge,  Stresses  In  ShellSi  Sprix^er-Verlag,  p,  580,  i960. 


TR-61-90 


a 


I 


:  1 
» »• 


;[ 

;[ 
;  I 


ii 

!  I 

:  r? 
!  ^ 

’[ 


I 

i 

1 

1 


DISTRIBUTION 

No.  Cvs 

HEADQUARTERS  USAF 
1  Hq  liSAF  lAFDDC),  XVash  Z5.  DC 

1  Hq  USAF  (AFOOA),  Wash  25.  DC 

1  Hq  USAF  (AFOlEi,  Wash  25,  D,; 

I  Hq  USAF  (AFRDR),  Wash  25,  DC 

1  Hq  USAF  UFCIN).  Wash  25,  DC 

I  USAF  Dep  IG  £or  Insp  (AFCDi-3-3).  Norton  AFD.  Calif 

1  Hq  USAF  (AFOSR.  Tech  Library),  Uldg  T-D.  Wash  25,  DC 

MAJOR  AIR  COMMANDS 
1  AFSC  (SCR).  Andrews  AFB,  Wash  25.  DC 

1  SAC  (DINC),  Oftutt  AFB.  Nebr 

1  TAC  (TDA).  Langley  AFB.  Va 

1  AFLC(MCW7),  Wright -Patterson  AFB,  Ohio 

1  ADC  (Ops  Anlys).  Ent  AFB,  Colorado  Springs,  Colo 

I  AUL,  Maxwell  AFB,  Ala 

1  USAFA,  United  States  Air  Force  Academy,  Colo 

AFSC  ORGANIZATIONS 

1  BSD  (BSR-VP-3),  AF  Unit  Post  Office,  l.os  Angeles  45,  Calif 
ASD,  Wright-Patterson  AFB,  Ohio 

2  (ASAPRL) 

1  (ASMCC) 

KIRTLAND  AFB  ORGANIZATIONS 
AFSWC,  Kirtland  AFB,  W  Mex 
1  (SWEH) 

26  (S<VOi) 

3  (SWRS) 

OTHER  AIR  FORCE  AGENCIES 

I  Director  USAF  Project  RAND,  ATTN:  RAND  Nuclear  Tnergy 

Division,  via  Air  Force  Liaison  Office.  The  RAND  Corporation, 
1700  Main  Street.  Santa  Monica.  Calif 

ARMY  ACTIVTIES 

»  Chief  of  Research  and  Develorm«-nt,  Department  of  the  Army, 

ATTN:  Special  Wcaponr  and  Air  Defense  Divistot,.  Wash  25,  DC 


I 


f 

\ 


1 


TR-61-90 


DISTRIBUTION  (con't) 

1  Director,  Ballistic  Research  Laboratories.  ATTN;  Librar>’. 

Aberdeen  Proving  Ground,  Md 

1  Commanding  Officer,  U3  Army  Engineers,  Research  and  Development 

Laboratories,  Ft.  Belvoir,  Va 

1  Office  of  the  Chief,  Corps  of  Engi:»e-*rs,  US  Army.  ATTN:  Pro¬ 

tective  Construction  Branch,  Wash  25,  DC 

NAVV^  ACTIVITIES 

1  Chief,  Bureau  of  Yards  and  Docks,  Departntent  of  the  Navy,  Wash  25, 

DC 

1  Commanding  Officer  and  Director,  Naval  Civil  Engineering  Labora¬ 

tory,  Port  Hueneme,  Calif 

OTHER  DOD  ACTIVITIES 

1  Chief,  Defense  Atomic  Support  Agency,  ATTN;  Document  l.lbrary 

Branch,  Wash  2b,  DC 

1  Commander,  Field  Command,  Defense  Atomic  Support  Agency, 

ATTN:  FCAG3,  Sandia  Base.  N  Mex 

10  ASTIA  iTIPDR),  Arlington  Hall  Sta,  Arlington  12.  Va 

1  Director,  Advanced  Research  Projects  Agency,  Department  of 

Defense,  The  Pentagon,  Wash  25,  DC 

*  OTHER 

1  Official  Record  Copy  (SWRS) 


,  2^1  as 

I  n  <4  M 

O  41  « 


h  q  o 

SS«>*=  4 

H  Z  M  U  ^ 


M  4  ^ 


as|  »s 

»lii|' 

rsspt 

A*4  4*  H  o  :; 

\  3  K  5 ».  J 
33"^ “t 

is^iis 

H  4  a  h  fii« 
3  S  o  •>  M  » 

h?s«f 

«  v4  U  O  VH  2 

l«i|; 

null 


u  5S*>S  ^  e  I> 

I  Sl|2  I 

SSolljgls  Ipibg 


£iiis§fi§S2saiss 


Mf4f^a\S  '^tm  IH  »H 


fiog 


X  •>.4  4  A  4  <14  U  3  M 

b,  o  M*>at««4i4  a  4>  9  V 

p  l"°:l|S8|  ,sifi 


"*§1  ®  ^aS“'So3i*"*.a°'g'S 

04*  4>  «»&.^l>  C^fi  M 

m  ^^!!P 

«3§e  fi|o^b'"3®*|M^5gS 

S«  •§  &  4>  «  4°  "§  >  *4  ^  5  XiJ  !3 

m. 

§2lf*  g3Ea8lH5sSI!!l 


I  I  |||f  I  |c|j| 

A  Sfi*^**  •  **  •*  •*  A  ^  S  •*  & 

i  IIIIIIIIillBIsiilfi 


V  V  >4 

si  nfIS 

a* 

J!5=ll 


mum 

3  >  Msc  o  &vt  u  e 
jin<  R.uua 


SxIsS  saE5 


-a.  3  3  8  .« 

0^>V  J«>-i5Sv4« 

'sPjlilfh 

S53''3°K1^S 

3a3fs|!f«|:i 

;p?p?3Ev 

8!;tR;aSa33 


